The FSI contribution to the observed B s decays into K + K and ty + K 
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Because at the tree level B s — > K + K~ is Cabibbo-triple suppressed, so its branching ratio should be smaller 
than that of B s — > -k + K~ . The measurements present a reversed ratio as R = B(B S — s- tt + K~)/B(B B — > 
K + K~) ~ 4.9/33. Therefore, It has been suggested that the transition B 3 — > K + K~ is dominated by the 
penguin mechanism, which is proportional to V c t,V* s . In this work, we show that an extra contribution from 
the final state interaction (FSI) to B s — s> K + K~ via sequential processes B s — > D^Di*^ — s> K + K~ is also 
£SJ , substantial and should be superposed on the penguin contribution. Indeed, taking into account of the FSI effects, 

the theoretical prediction on R is well consistent with the data. 
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I. INTRODUCTION 



_ It is generally believed that the hierarchy of the Cabibbo-Kabayashi-Maskawa (CKM) matrix elements determines the mag- 
nitudes of weak decays, and off-diagonal matrix elements would bring up an order suppression. For example, the ratio 
of B(B S -> D+tt')/B(B s -> DfK T ) = 3.2 x 10~ 3 /3.0 x 10~ 4 has been measured [1, 2], and it is roughly deter- 
mined by the ratio of the CKM matrix elements V u d/V us . Thus one usually categorizes the weak decays according to their 
CKM structures as the Cabibbo-favored; Cabibbo-suppressed and even the Cabibbo-double or triple-suppressed. However, 
Q-[ the newly observed modes B s K + K~ and B s — > ir + K~ obviously do not follow the rule, namely at first look, the 
branching ratio of B s — > K + K~ should be smaller than that of B s — > ir + K~ by |Kis/Kid| 2 > by contraries, the datum is 
K = B(B S -> n+K-)/B(B s -> K+K~) = 4.9 x 10 -6 / 3 - 3 x 1CT 5 [3-5]. If only the tree diagrams are taken into account, 
this reversion would compose an "anomaly". Descotes-Genon et al. [6] carefully analyzed the transitions of B s — > K°K° 
and K + K~ through flavor symmetries and QCD factorization, and pointed a potential conflict between the QCD prediction on 
^ ■ B s K + K~ and data. By analyzing the transition mechanism, Cheng and Chua [7] determined that the main contribution 

■ to B s — > K + K~ comes from the penguin diagram, which is proportional to V c bV* s . Ah et al. also calculated the branching 
^ ! ratios of B s -> tt+K~ and B s -> K + K~ in terms of the Perturbative QCD (PQCD) to LO and NLO [8, 9] and the authors of 

■ Ref. [10] did calculations in SCET. Their results roughly were consistent with the data available then, so the "conflict" seemed 
• | resolved. Even though the theoretical uncertainties in all the calculations are not fully controlled, it is noted that the obtained 

central values are not sufficiently large to make up the data. It implies that there must be some mechanisms to remarkably 
' enhance the branching ratio of B s — > K + K~ . 
t—( . Looking at the central values they obtained and the resultant ratio of B(B S — > it + K~)/B(B s — > K + K~), one can find that 
I* 1 the calculated B(B S — > ir + K~) is close to the data, however, the central values of B(B S — > K + K~) calculated in various 
approaches are smaller than the newly measured data [4, 5]. 

Based on this observation, we suggest that the final state interaction (FSI) in B s decays may greatly enhance the branching 
ratio of B s —> K + K~ but not much for B s — > ir + K~ . In fact, in the energy regions of 6-quark and c-quark most such anomalies 
can be naturally explained by considering the role of FSI. For example, a simple quark diagram-analysis tells that the branching 
ratio of D° — > K°K° is almost zero, but its measured value is comparable with that of D° — > K + K~, which is large. This can 

be understood by considering the sequential process D° — > K + K~ — > K°K° and the later step is a hadronic scattering [11]. 

(*}+ (*}— 

The Particle Data Group (PDG) [3] tells us that the D s D s are the dominant hadronic decay modes of B s , therefore it 
implies that the sequential decays B s Dj*' + Di"' —> K + K~ would compose an important contribution to the observed 
B s — > K + K~ . The first step of B s — > Ds + is only suppressed by Vb c and the process does not suffer form a color- 

suppression. Moreover it is also isospin-conserved mode, even though the weak interaction does not require isospin conservation, 
it still may be more favorable than the isospin violated ones. Thus, one can understand why B s — > Ds Dl is dominant. 
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Then let us look at the next step. The FSI is a hadronic scattering process, which is completely governed by strong interaction, so 
that the isospin must be conserved. The isospin of D s is zero, thus, the final state of the scattering is required to be zero. The 
isospin of A'-mesons is 1/2, thus the KK -states can be either isospin or 1, therefore the inelastic scattering D\ ' D\ ' — > 
K + K~ is allowed. By contraries, isospin of pion is 1, thus, the system of ir + K~ does not contain an isospin component, 
thus the scattering D\ D s — > n + K~ is forbidden. Definitely, one can expect a substantial contribution from B s — > 
£>(*)+£>(*)- K+K-, which can much enhance the branching ratio of B s — > K + K in comparison with B s — > tt + K . 

It is worth pointing out that on the other hand, the decay B s — > ir + K~ also receives a contribution from the FSI via B s — > 
£)(*)+£)(*' _^ tt + K ~ . However, the first step process B s —> £)(*)+ Ds*^~ is Cabibbo-double suppressed by V c bV* d , so is 

( S|e ^ -j- f S|i ") _|_ 

much less than B s — > D s ' D s , therefore the FSI does not contribute much to the branching ratio of B s — > n + K~. 

Notice that, the direct process B s — > K + K~ via penguin diagrams and the sequential process with FSI B s —t — > 
K + K~ have the same initial and final states, moreover their amplitudes are of the same order of magnitude, so the two con- 
tributions interfere. In fact, their contributions are not directly experimentally, but theoretically distinguishable. Moreover, the 
strong scattering would have a real and an imaginary parts (see below, for the calculations of the triangle diagrams), thus a phase 
is resulted. In most calculations of the FSI effects, only the absorptive (imaginary) part is kept, the reason is that one may argue 
that the absorptive part might be dominant or at most the dispersive and absorptive parts have close magnitudes. In that case, the 
absorptive part is imaginary while the penguin contribution is real, so that the two contributions can be added up at the rate level 
(amplitude square). Our final results indicate that while taking into account the new contribution, the theoretical predictions are 
indeed close to the new data. 

In this work, we calculate the amplitudes of the decay channels of B s — > Z?i*' 1 )Z?i*' 1 and the direct decay channels 
of B s — > K + (ir + )K~ by the factorization approach [12, 13]. Then we use the effective SU(4) Lagrangian [14, 15] to 
determine the vector-pseudoscalar-pseudoscalar and vector-vector-pseudoscalar vertices for calculating the FSI amplitude of 
£>W (£>(*) )£W -> K+{tt+)K-. The paper is organized as follows, after the introduction, we formulate the weak decays and 
the re-scattering processes in section II, our numerical results are presented along with all necessary inputs in section III. The 
last section is devoted to discussion and conclusion. 



Because the measurements on the decay widths of B s — s- D s * D s * '~ and B s -> D^ + D s (B s -> D^~ D 8 *' + ) are 
not accurate yet, as for most of the above channels there are only upper bounds, so that we are going to directly calculate the 
transition amplitudes based on the quark diagrams. Even though this strategy might bring up certain theoretical uncertainties, it 
does not break our qualitative conclusion at all. 

For calculating the transition amplitudes of B s — > Ds*^ + )D S *^ , one needs to employ the effective hamiltonian at 

the quark level. With the operator product expansion (OPE), the effective Hamiltonian was explicitly presented in Ref. [16]. 
At the tree level, from the effective Lagrangian one can notice that B s — > K + K~ is suppressed by V u bV* s i.e triple-Cabibbo 
suppressed, comparing with B s — > tt + K~, which is double-Cabibbo suppressed by V u b- Therefore Cheng and Chua decided 
that the direct transition B s — > K + K~ is obviously dominated by the penguin diagram whose CKM structure is V c bV* s . The 
transition B s — s- D s *^ + Z)W _ is also double-Cabibbo suppressed. Instead, B s D)* Ds is proportional to V c bV* s and it 
is a tree process. Now let us compare the sequential processes B s — s- D^Di,*^ — s- K + K~ with the penguin contribution. The 
two reactions are of the same CKM structure, but the penguin undergoes a loop suppression about a s /ir ~ 0.06, whereas the 
strong scattering Ds D s ' — > J2i where stands as any possible final states allowed by symmetry and energy-momentum 

conservation. K + K~ is only one of the possible channels, and its probability is proportional to (Ds*" 1 D s *^ \H e f f\K K), which 
is what we are going to calculate in this work. This is a suppression factor because the total probability to all channels is 1 . 
Therefor, roughly, we notice that the penguin is loop-suppressed and the sequential process is also suppressed by the probability, 
thus the two modes compete and may have a similar order of magnitude. Concretely, we need to calculate them. The explicit 
calculation on the penguin contribution can be found in [7-9], thus we will use their numbers and only consider the contribution 
from the sequential processes. 

Applying the effective hamiltonian at the quark level to the hadron states, the hadronic matrix elements can be parameterized 



II. THE THEORETICAL EVALUATIONS OF THE BRANCHING RATIOS 



A. Weak decays in factorization approach 



as [17]: 



(0\J,\P(Pi)) 
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{P{ P2 )\J tl \B s {p)) = P» ~ 2 " F 1 (q*)+ -»° , - q^Fotf), 



<1 



m. 



<1 



(V(p 2 ,e)\J^\B s (p)) = lJ ~r \i^ a pP a q fi A v (q 2 ) + ( mBs +m v ) 2 g^A 1 (q 2 ) ^_A 2 (g 2 ) 

m Bs + m v I m Bs + m v 

-2m v (m Bs + m v )^[A 3 (q 2 ) - A (q 2 )}}, (1) 

where J M = gi7 M (l - Js)Q2, = (Pi + V2)p. and q^ = (p x - p 2 ) t _ l . 
With Eq. (1), we write down the amplitudes of B s — > D s D s : 

A{B s (p) -+ Dt( Pl )D-{p 2 )) = -^V cb V c * sai f Ds (rn 2 Bs - m 2 D jF* ° D » (p 2 ), (2a) 
A(B S ( P ) -> D s *+( Pl )Dr(P2)) = ^VctV'oi/o.mo; mfl ^ mg . {*w(pi + 2pa) a pf (rf) + (m Bs 

+m D ,y 9 » v A B 1 ° D '°(j>{) ( Pl + 2p 2 ),( Pl + 2p 2 )X sD: (pi) - {pi + 2p2) ^ \(m Bs 



+m_ D ») 2 J 4f s£ ' s (p 2 ) - (m| s - m%,)A2 aDs (p 2 ) - 2m D *(m Bs + m D *)A^ sDs (p\) X 



(2b) 

and the amplitudes of _B S — » Z?(*' + _Ds*' ) ~ read as 



A(B s (p) -+ D+{ Pl )D-(p 2 )) = - l -±V cb V L * d aif D (m 2 Bs - m 2 Dt )F*' D ° (p 2 ), (2c) 

MB.{p) -> D^^D*- (pa)) = ^y cb y c >i/j>m . {»W(Pi + 2p 2 ) a p?< s ^ (p 2 ) + (m Ss 

+m D .)V^(P?) - (Pi +2p 2 ),(p! +2p 2 )^f ^*(p 2 ) _ (Pi+ 2 ^W r (mBs 

Pi L 

+TO r ,») 2 -4f 3£ ' s (Pi) - (m| a - m 2 D ,)A2 sDs (p 2 ) - 2m D *(m Ba + m D *)A^ sDs (p 2 ) |, 

(2d) 

the amplitude of 7r + A' - at the tree level is: 

A dilcct (B s (p) -> ^+(pi)A-(p 2 )) = -^V^KViACmi. - )A s ^(p 2 ), (2e) 

where ai is a proper combination of the Wilson coefficients in the effective hamiltonian [16, 18, 19]. Fo(q 2 ) and Ay,i, 2 {q 2 ) 
are the form factors to be detrmined. In this work, due to lack of accurate data, we use the form factors obtained by fitting 
the data of the decays of B meson. It is a reasonable approximation because the processes B s —> D* and B —> D* have the 
same topological structure, and the flavor SU (3) symmetry for light quarks (it, d and s) would lead to the same form factors, 
namely the difference between the form factors for different light-quark flavors would be proportional to an SU(3) breaking, 
which is small for the effective vertices as well known. At least such small difference would not overtake the errors caused by 
experimental measurements and theoretical uncertainties for evaluating the non-perturbtive effects. 

Moreover, a symmetry analysis indicates that the sequential process B s —> Z)* + D~ (£)+ -D* - ) — > K + K~ is forbidden by 
angular-momentum conservation. 

Taking the three-parameter form, the form factors are written as [17]: 



where a, b and F(Q) are the three parameters and their values are listed in Table I. 
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F F(0) a b 


F F(0) a b 


Fo D 0.67 0.65 0.00 
V BD " 0.75 1.29 0.45 
Af D * 0.63 0.65 0.02 


Fo n 0.25 0.84 0.10 
Aq D * 0.64 1.30 0.31 
Af D " 0.61 1.14 0.52 



TABLE I: The parameters given in Ref. [17]. 




(d) 



(e) 



(f) 



FIG. 1: The QCD one-particle exchange process in Z?s*' + —¥ K + K 



B. Evaluation of FSI effects 



Now let us turn to evaluate the long-distance effects at hadron level. For the effective vertices, the flavor- SU (4) symmetry is 
assumed. The coupling of pseudoscalar and vector meson is 



1. 



Co = Ti(d^d^) - -Tr(Ft F^), 



(4) 



where F^ v = d^V u — d v V^, and <j> and V represent the 4x4 pseudoscalar and vector meson matrices in SU (4), respectively: 
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For the gauge invariance, covariant derivatives replace the regular ones: 



(6) 



F, 



(7) 
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Now we are ready to write down the relevant terms in the pseudoscalar-vector coupling: 

C = £ + igTv{d»(t>[4>, V v ]) + g'e^d^d^V^ + 



(8) 



where the third term on the right-hand side of Eq. (8) is involved additionally as a vector- vector-pseudoscalar coupling [23]. In 
the process that Di*^ + Di*^~ re-scatter into K + K~, the corresponding Lagrangian is: 



£ 
C 

£>KD*D 



{ 



d: 



KD'D — igKD'D | — sfi 



>{ 

gKD»D*£ 



D 



D 



Sp 



}■ 



Y 



(9) 



whereas the effective vertices for irDD* and irD*D* can be found in Ref. [14]. 

In the Eq. (9), the values of the coupling constants should be obtained by fitting the experiment data. However, it is noticed 
that in previous literature those coupling constants are still not well fixed yet (see Table II), so in this work, we use an average 
value for each coupling constant, and the uncertainty is considered as a systematical error, which might be attributed to our input 
parameters. For example, there are no available data for determining the coupling constants gKD s D* and gKD*D, we are going 
to fix them based on the SU (4) symmetry, which tells us that: gKD s D* = gKD' d = VZgpTrrr = VSg^DD/^, and by fitting the 
data, g pn7T = 8.8 and g^DD = 7.9 are set (see Table II). We may have two different values for gx d s d* and gKD* d- 6.0 and 
9.4, but by our strategy we adopt an average value for gKD s D* and gKD*D as: gKD s D* = gKD*D = 7.7 ± 1.7. We use the 
same method to get the values for other relevant coupling constants as long as there are no data to directly fix them, and retain 
the errors. Then we have: 



gi<D s D" = gi<D" s D = V2g P inr — —^-g^DD = 7.7 ± 1.7, 
gKD*D* = ~^^' D * D = 3.2 ± 0.8GcV _1 , 
g-xD*D> = 9.0 ± 0.1 GcV-\ g vDD . = 4.4. 



(10) 



To evaluate the FSI effects, we calculate the absorptive part (see above arguments) of the one-particle-exchange triangle 
diagrams. Fig. 1 shows the diagrams of B s —> d{*^ + d{*^ —> K + K~ by exchanging D°(*\ For calculating the absorptive 
part of the triangle, with the Cutkosky cutting rule [20, 21], the related amplitudes are: 



B s (p) Dt( Pl )D:{ P2 ) -+ K+(p 3 )K-( P4 ) 



dp 1 dp 2 (2n) 4 S(p-p 1 -p 2 )A 



D t(?l) D s (P2) {-gKD„D») ~i(pi +P3)fi (-gKD 3 D>) ~ i(p2 + Pa)u 



g^v 



i D o, 



-2 2 - 7r2 ( l ? 2 > m i>o»), 



(11a) 



A 2 B s (p) D*( Pl )D*(p 2 ) -> K + ( P3 )K-( Pi ) 



- I dp 1 dp 2 (2Tr) i 5(p-p 1 -p 2 )A 



9pi 



PinPit 



m 



D" 



g V a 



B s 

P 2v P 2cr 



m 



D • 



Dt*{ Pl )D-*( P2 ) ( 9kd*d) 
i 



KQ+P3)c ( 9kd*d) 



m 



D" 



(lib) 



A; 



B s (p) -> D*+( Pl )Dr( P2 ) -> K+(p 3 )K-( Pi ) 



- / dp^p^ir) S(p- pi - p 2 )A 



B„ 



Df*( Pl )D-*(p 2 ) {ig K D;D*)Hpi*)(-iq a )*™ aP (igKD;D*)HP2e)(iq K ) 



9Pp 



<- D o* 



gpu) 



PlpPlu 



D ' 



g v \ 



P 2u P 2\ 



iir. 



-F 2 {q 2 ,m 2 D , ) 



(11c) 
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instead, for the decay channel B s —> D^ + Di*^ —> tt + K , we have: 



B s {p) -+ D + ( Pl )Dj(p 2 ) -> 7r + (p 3 )tf-(p 4 ) 



dp 1 dp 2 (2TT) 4 S( P - Pl - P2 )A B s ^ D + { Vl )D-(p 2 ) (-g« DD .) 



i(pi +P3)^ 



9^ 



q»q» 



i< D o, 



-= g • ?72 (9 2 ) m IJO«)) 



(-9kd s d-) -i(P2+P4)i 

(Hd) 



-As 



B a (p) -> D* + { Pl )D*-(p 2 ) -> 7r+(p 3 )^-(P4) 



= ^ / dpidp 2 (2-K) A 5(p- pi -p 2 )A B s -> D + *(p 1 )D s *{p 2 ) (-g^DD') -i(q+P3)i (-9kd*d) -i{q-pi) c 



9rf 



g 2 - m^o 



J" 2 (g 2 ,m|,o), 



(He) 



-4a 



B fl (p) -> D* + (pi)D;-(p2) 7r+(p 3 )A'-(p4) 



^ / dp 1 dp 2 (2 7 r) 4 (5(p-p 1 -p 2 M B s ^ J D+*(p 1 )£) s -*(p 2 ) (iff 7rD ./3.)(-ipi CT )(-ig»)e CTWa ^(iffx D jzj«)(-^ 2 {)(^ K ) 



Xe 



9Pp 



qpq P 



9tioj 



9v\ 



V2vV2\ 



■111-, 



q 2 - m 2 D0 , 



T 2 {q 2 ,m 2 D , ), 



(Hf) 



where dp = dp 3 / ((2ir) 3 2E), q is the momentum of the exchanged D°(*) meson and a dipole form factor J-(q 2 , m 2 ) — (A 2 
m 2 ) 2 /(A 2 — q 2 ) 2 with A = m + (3Kqcd 1 is introduced to compensate its off-shell effect [22]. 



g value 


g value 


g p ^ 8.8 [14] 


gDDii) 


7.9 [22] 


g D * m , 4.2GeV~ 1 [22], 7.7 [14] 


gwD'D' 


8.9GeV _1 [22], 9.1 [24] 


g lDD . 4.4 [23] 







TABLE II: The values of the coupling constants involved in our work. 



III. NUMERICAL RESULT 



With the amplitudes given in Eq. (1 1), we can easily calculate the decay width of the sequential processes B s — > Dg Dg — > 
K+K~. We take the relevant parameters from PDG [3] as: m Bs = 5.366 GeV, m Ds = 1.968 GcV, m D * = 2.112 GeV, 
m D o — 1.864 GcV, m D o* — 2.007 GeV, m D ± = 1.869 GcV, m D *± = 2.01 GcV; and the value of the weak decay constants, 
such as fir Jk fn Id, can be found in Ref. [25]; ai = 1.14 is set in our numerical computations [13]. The total amplitudes are 



A(B S -> K + K~) 



A dilcct (B s -> K+K-) + A\ (B s -> D^D- K + K~) + A 2 {B S -> 
MsCB, -> D*+D*- -> if+ii"-) = 4.1±i;S x 10~ 8 , 



K+K~) 



(12) 



^dircct^ ^ tt+A'-) +^(5, -> D + D~ -> TT+iT") + .4 2 (B S -> 
+-4 3 (£s -> D* + D* s - -> Tr+if-) = 29.5±g:i x 10~ 9 , 



(13) 



Here, m denotes the mass of the exchanged meson and /3 is a phenomenological parameter. Aq<j £> is taken as 220 MeV. 
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where A ducct (B s — > K + K ) is the penguin contribution [7]. So, the branching ratios are: 



B(B a ->■ K+K-) 
B(B S tt+A' - ) 



1 B S ^K+K- 
Ttot 

Fb s ^ttK 



1 



327r 2 r tot m 2 B 



A to t(B s -> KK) dtt = 13.6+7 5 7 5 x 10~ 6 , 



1 



PK 



327r 2 r to t m% 



Aot(S s . -> At^ 



rfO = 7.3i°; 1 x 10"°, 



(14) 



where the errors are systematical, originating from the uncertainty of the coupling constants. 



TABLE III: Various approaches predictions on B(B S —>■ K + K ) (in units of 10 6 ). 
and the ratio TZ = B(B S -> ty+K-)/B(B s -> K+K~) 





QCDF [7] 


pQCD (LO) [8] 


pQCD (NLO) [9] 


SCET [10] 


FSI 


FSI+pQCD (NLO) 


Experiment [4, 5] 


B 


r,r r, + 12.7+12.5 

— 7.2—9.1 


13.6± 8 5 i 


15.6i!;J 


18. 2 ±6. 7 ±1.1 ±0.5 


13.6± 7 5 7 5 


29.2+g 5 7 8 


33 ±9 


TZ" 


0.194 


0.360 


0.314 


0.269 


0.360 


0.167 


0.148 



"Since the main contribution for the transition B s — > n+K comes from the tree diagram and all predictions made in various models on this channel are 
close to each other, we use the data B(B S — > tt + K~) = 4.9 X 10 — 6 as input in our calculations. 



In Table III we list the theoretical predictions on B(B S — > K + K~) in various approaches. It is noted that most of the predicted 
cental values are lower than the newly measured value [7-10] as long as the contribution from FSI is not included. We add the 
contributions from FSI to that calculated in pQCD (NLO), then one can find that the resultant central value is consistent with 
data [30]: B(B S -> K+K~) = (3.3 ± 0.9) x 10~ 5 within la. 

The deviations of our theoretical prediction from the data might come from the loophole in our calculation. As indicated above 
we only consider the absorptive part of the hadronic triangle. Indeed the dispersive part may also make substantial contributions 
[26]. As suggested in literature, the dispersive contribution should be smaller than the absorptive one (it is consistent with the 
general principle of the quantum field theory), or at most has the same magnitude as that of the absorptive part. If the contribution 
of the dispersive part is indeed of the same order as the absorptive part, then taking it into account, we may have a result, which 
is even closer to the data. 

In general, even though we cannot precisely re-produce the experimental data, we can confirm ourselves that the FSI is 
important and non-negligible for understanding the B(B S — > K + K~) > B(B S — > ir + K~) "conflict". 



IV. CONCLUSION AND DISCUSSION 



As expected, the LHCb is extensively aiming on study of the B-physics, especially to look for some "anomalies" in experi- 
ments, which need high statistics and precise measurements. For i? s 's charmless non-leptonic two-body decay, the early MC 
studies show that, nearly 37K B s — > K + K~ signals will be seen at 2 fb^ 1 integrated luminosity [27]. On the other hand, 
314 ± 27 B s —> ir + K~ signals were observed based on the 2010 data with its integrated luminosity 0.35 fb^ 1 [28]. We believe 
that the statistics of the B s decay into K + K~ and 7r + /\ ~ is sufficient to draw a definite decision about their branching ratios. 
The latest result reported by the LHCb collaboration shows that the branching ratios have been measured as [29]: 

B(H b -> F) = fw b N(H b ->■ F) e ICC (H' b -» F') e PID (F') 

B{H' b -> F>) f Hb N{H' b -> F') e rcc (H b -> F) e PID (F) ' 1 

where the /„o is the possibility of b quark hadronizing into hadron H, N is the observed number of signals for certain decay 

modes, e rcc is the efficiency of the reconstruction excluding the particle identification (PID) cuts and epio is just the efficiency 
of PID cuts. 

It is noted that in Table III, the experimental data are taken from Refs. [4, 5], but the 2012 data of PDG indicate that the 
branching ratio of B s -> K + K~ is (2.64 ± 0.28) x 10~ 5 [30], which is smaller than the data in Ref. [4, 5]. Moreover, the 
LHCb collaboration reports that with the 0.37/6 -1 2011 data, the branching ratios of B s — > K + K~ and B s — > tt + K~ are 
experimentally determined as: 



B(B S -> K+K-) = (23.0 ± 0.7 ± 2.3) x 10~ 6 , 
B(B S -> 7T+K-) = (5.4 ± 0.4 ± 0.6) x 10" 6 , 



(16) 
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where the former uncertainties are statistical and the later one is systematical, which include the uncertainties of PID calibration, 
final state radiation with soft gamma, signal shape used for fitting, and the impact of the background: the additional three-body 
background 2 , the combinatorial background and the cross-feed background 3 . Since B s -> K+K- and B s D^D^ -> 
K + K~ have the same final states, it means we cannot single out the contributions of FSI by simply measuring the cross section 
in experiments. Even though our theoretical prediction on B(B S K + K~) presented in Table III is slightly above the value of 
the LHCb measurements, it is still consistent with the data within the experimental error tolerance. Therefore we are expecting 
more precise measurements in the future. 

Cheng and Chua suggested that the decay B s — > K + K~ is dominated by the penguin diagram, which does not suffer from 
large CKM suppression [7]. In that scenario the "conflict" B(B S — > K + K~) > B(B S —> ir + K~) could be partly explained. 

In our work, we show that the FSI definitely makes an important contribution because B s Dl ' D\ ' are dominant decay 

(^"l _|_ _ _|_ 

modes of B s as confirmed by the data and the scattering D), D s — > K + K~ is allowed by all the symmetry requirements. 
Therefore, we might consider that the penguin contribution and the FSI effects are in parallel to contribute to the decay B s — > 
K + K~ . A Actually, there are several phenomenological parameters in the model for evaluating the FSI effects, which were 
obtained by fitting the earlier data of heavy flavor processes. One thing can be sure that the FSI should exist and contribute to the 
process B s — > K + K~~ , but how significant it would be is determined by both requirement of experimental measurements and 
theoretical estimate. If the data on B(B S — > K + K~ ) are indeed going down, the FSI effects for B s — > K + K~ would be less 
significant, by that situation, one can further restrict the involved phenomenological parameters at this energy region. On other 
aspect, the errors of the measurements are still too large to make a definite conclusion yet, so we are expecting more precise 
measurements not only from LHCb, but also the future super-B factory. 

Definitely the FSI effects also play roles in other similar decay channels, such as B ± — > K^ui. Therefore careful studies 
on such modes with the same theoretical framework would be helpful for more accurately evaluating the FSI effects. The PDG 
indicates that, B(B ± -> A^w) /B(B ± -> tt ± uj) = (6.7 x l(T 6 )/(6.9 x 1CT 6 ), so one would expect that B(B -> Ku) is also 
enhanced by the FSI. Moreover, a measurement on the polarization of lu might be helpful to gain more information about the 
role of the FSI mechanism, and it will be studied in our coming work. 

From the experimental aspect, as the FSI is a scattering fully governed by the strong interaction, its proper time is too short 
to be measured in the LHCb detector. Although the LHCb can well reconstruct the events B s — > K + (tt + )K~ and identify 
the it + K~ purely with the PID cuts [31], it is impossible to distinguish between the direct decays B s — > K + (tt + )K~ and 
the sequential ones. Considering the good ability of LHCb for tracking charged particles, we believe that, B ± — > A'(7r) ± o; is 
another good channel to study FSI. The decay rates of B(B — > D°D^*' + ) have been well measured as ~ 1.8 ± 0.2 x 10 -2 , 

and we can make Mont-Carlo simulations on the re-scattering of D°Dg + into K(ir) + uj to complete the theoretical scenario. 
The abnormal ratio of the two channels B(B ± -> K ± uj)/B(B ± -> tt^w) = (6.7 x l(T 6 )/(6.9 x 10" 6 ) [30] can also be 
explained as the FSI contribution as well 5 . With B ± — > £>A'(7r) ± measurement [32] and considering the efficiency, about 0.8K 
B + —> K(it) + uj can been seen in the 2011 database and it implies that the FSI contribution to the polarization of the vector 
meson ui might be distinguished from that of the direct decay. Comparing more accurate data with our theoretical calculation, 
one can expect to gain more knowledge on the FSI effects, for example how to determine the parameters in the dipole form 
factor etc. 

One object can be recommended for pining down the role of the FSI effects. As well known that the direct CP violation is 
proportional to sin(o:i — a^) ■ sin(</>i — <f>%), where the subscripts 1 and 2 refer to two different routes and a, <fi are the weak and 
strong phases respectively. Therefore, if a direct CP violation exists in B s decays, there must at least be two different routes to 
the final states which possess different weak and strong phases. In fact, the penguin and tree diagrams have different strong and 
weak phases. The weak phase of the tree diagram is coming from V u bV* s and does not posses a strong phase. Instead, the weak 
phase of the penguin is from V c bV* s and the strong phase is due to existence of the absorptive part of the loop. In our scenario 

the sequential processes B S (B S ) — > DgDi — > K + K~ which have the same weak phase V c bV* s , but different strong phase 
from the penguin. Because of the extra contribution, the amplitude would become 

]\J — J\/jtree e ia! _|_ j^j-penguin e i(a 2 +4>l) _|_ ]\/[seq e i(a 2 +4>2) 

where the last term is the FSI contribution. 

Thus the interference between the tree diagram with a sum of the penguin and sequential processes would be different from 
the situation where the tree diagram only interferes with the penguin. Therefore by measuring the CP violation of B(B S — > 



2 Miss or misidentify the pion or kaon in final state. 

3 The uncertainties from the distribution of signal in data and simulation. 

4 In this work we do not consider the penguin contribution, but only that of the FSI effects to the decay. We acknowledge the possible contribution from the 
penguin diagram as well, further study involving both of them and moreover their interference will be made in our later works. 

5 B(B -»• K (7r)w) = 6.7(6.9) x 10" 6 , at the same time, the Cabibbo suppression determines that the branching ratio of Kcu final state should be smaller. 
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K + K~) — B(B S — > K + K~~), one can distinguish between the penguin contributions and the FSI effects. But since the 
experimental errors and theoretical uncertainties are not well controlled so far, we cannot make a more definite prediction on the 
CP violation yet, we would wait for more precise data to be available and then continue our theoretical computations. 

Moreover, according to the above arguments, there does not exist a tree diagram for B s — > K°K°. Consequently, the direct 
transition B s — > K°K° is uniquely determined by the penguin diagram and FSI effects. Because the tree contribution to 
B s — > K + K~ is very small, the rate of B s — > K°K° must be close to that of B s — > K + K~. However the interference 
between penguin (neglecting the Cabibbo-suppressed penguin diagrams where i-quark and u-quark are intermediate agents) 
and FSI does not lead to a direct CP violation at all, because they have the same weak phase. This can be confirmed by future 
experiments. 
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